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ABSTRACT 

4U 1626—67 is an accretion powered X-ray pulsar that shows remarkably stable 
X-ray luminosity above hours timescale and gradual intensity variation on a few years 
timescale. Unlike the other high magnetic field binary X-ray pulsars, the spin-up or 
spin-down rate of the neutron star is also very stable upto several years. Here we report 
a significant increase in the X-ray intensity in the long term RXTE-ASM light curve 
of 4U 1626—67. Similar enhancement in the X-ray flux has also been detected in the 
Swift-JiAT light curve. The increase in the X-ray flux took place over a long period of 
about 100 days and there appears to be two episodes of flux enhancement. We have 
investigated the spectral and timing features of 4U 1626—67 during its current state of 
enhanced flux emission with data obtained from the Proportional Counter Array and 
the High-Energy X-ray Timing Explorer on board the Rossi X-ray Timing Explorer. 
We report the detection of a torque reversal to spin-up in 4U 1626-67. The source has 
entered a new spin-up phase with a spin-up rate of 4.02(5) x 10"^'^ Hz s~^. The present 
spin-up rate is almost half of the earlier spin-up and spin-down trends. A significant 
excess in soft X-ray photon emission is observed during the enhanced flux state, which 
is similar to the energy spectrum obtained during the spin-up era of the pulsar before 
1990. 4U 1626—67 is a unique accretion powered X-ray pulsar in which quasi periodic 
oscillations have been consistently observed over the past ^^20 years. But during the 
recent observations, we did not detect a QPO at frequencies as observed in earlier 
observations. Instead, we report detection of a significant broadening in the wings of 
the 130 mHz peak and a change in the shape of the continuum of the power spectrum. 
These results indicate that the flux enhancement is not a simple case of increased mass 
accretion rate, but there is also a change in the accretion geometry in the vicinity of 
the neutron star. 
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1 INTRODUCTION 

4U 1626—67 is an accretion powered X-ray pulsar which 
was discovered with the Uhuru satellite in 1972 (Giacconi 
et al. 1972). The 7.7s X-ray pulsations were discovered by 
Rappaport et al. (1977) during SAS3 observations and have 
been extensively monitored since then. Initially the pulsar 
was found to be spinning up with a characteristic timescale 
of ^5000 years. It underwent an abrupt torque reversal 
in 1990 and has been steadily spinning down since then 
(Chakrabarty et al. 1997, Krauss et al. 2007). Optical coun- 
terpart of the pulsar was identified as KZ TrA, a faint blue 
star with little reddening (McClintock et al. 1977; Bradt & 
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McClintock 1983). Pulsations at 130 mHz have also been de- 
tected from the reprocessed optical emission (Middleditch et 
al. 1981). Multiple sidebands separated by about 0.4 mHz, 
are present in the power spectrum of the optical light curve, 
which was interpreted as signature of an orbital period of 
~42 min (Middleditch et al. 1981; Chakrabarty et al. 2001). 
But despite extensive searches, an orbital motion of the neu- 
tron star (usually observed in binary pulsars in the form of 
Doppler shift in the pulse period, or delay in the pulse arrival 
time) has never been detected in this source (Rappaport et 
al. 1977; Levine et al. 1988; Jain et al. 2007). It is believed 
that the neutron star most probably, has a low mass binary 
companion and is in a nearly face-on orbit. 

The X-ray spectrum of 4U 1626—67 was first studied in 
detail with the ASCA spectrometers (Angelini et al. 1995) 
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and later with instruments on board Beppo-SAX (Owens 
et al. 1997; Orlandini et al. 1998), Chandra (Schulz et al. 
2001) and XMM- Newton (Krauss et al. 2007). The X-ray 
continuum spectrum is evidently correlated with the torque 
state. During the 1977—1990 spin- up phase, the spectrum 
was well described by an absorbed blackbody, a power law 
and a high energy cutoff (Pravdo et al. 1979; Kii et al. 1986). 
But after the torque reversal in 1990, the time-averaged X- 
ray spectrum was found to be relatively harder (Vaughan & 
Kitamoto 1997, Yi & Vislmiac 1999). Orlandini ct al. (1998) 
reported an absorption feature at ~37 keV which they at- 
tributed to be the cyclotron absorption line. The spectrum 
is rich in bright hydrogcn-likc and hcliurn-likc oxygen and 
neon emission lines. The blue and red Doppler shifted emis- 
sion lines in the X-ray spectrum indicate that the lines orig- 
inate in the accretion disk. However, presence of the red and 
blue Doppler shifted emission lines from the accretion disk 
is difficult to explain if the binary system is face-on. Pres- 
ence of a metal rich accretion disk and absence of orbital 
motion led Jain et al. (2007) to suggest that 4U 1626—67 
is a candidate for a neutron star with a supernova fallback 
accretion disk. 

An important signature of the presence of an accretion 
disk is the detection of Quasi Periodic Oscillations (QPOs) 
in the power spectrum of the X-ray pulsar. QPOs are be- 
lieved to arise due to inhomogeneities in the inner accretion 
disk. This aperiodic X-ray variability has been detected in 
several X-ray sources, such as, Cen X-3 (Takeshima et al. 
1991; Raichur & Paul 2008), XTE J1858-h034 (Paul & Rao 
1998; Mukherjee et al. 2006), EXO 2030-^375 (Angelini, 
Stella & Parmar 1989), XTE JOlll.2-7317 (Kaur et al. 
2007) and 3A 0535-^262 (Finger, Wilson & Harmon 1996). 
4U 1626—67 is a unique X-ray pulsar in which QPOs have 
been consistently detected in the X-ray and optical observa- 
tions. The QPO frequency has been evolving over the past 
22 years (Kaur et al. 2008). The QPO central frequency was 
found to be increasing during the spin-up phase of the pul- 
sar, whereas, a negative trend was observed after the torque 
reversal. During the spin-up phase of the source, QPOs were 
detected at a frequency of 36 mHz with the EXOSAT data 
(Kaur ct al. 2008) and at 40 mHz from the Ginga observa- 
tions (Shinoda et al. 1990). During the spin-down era of the 
source, QPOs have been reported to occur at a higher fre- 
quency of about 48 mHz from ASCA, Beppo-SAX, RXTE, 
and XMM-Newton (Angelini et al. 1995; Owens et al. 1997; 
Kommers ct al. 1998; Krauss et al. 2007). Large amplitude 
0.3—1.2 mHz QPOs were also detected in reprocessed op- 
tical emission from both ground based and Hubble Space 
Telescope observations (Chakrabarty et al. 2001). 

Recent observations with the monitoring X-ray instru- 
ments on board the RXTE and Swift satellites, show an in- 
crease in the X-ray intensity of 4U 1626—67. In the present 
work, wc have investigated the spectral and temporal prop- 
erties of 4U 1626—67 during the increased flux state. In par- 
ticular, we have searched for pulsations and the QPO feature 
in the power spectrum of 4U 1626—67 which was persistently 
present in the earlier observations. We have also studied the 
energy spectrum of the source during the current phase and 
have compared it with the earlier known spectral features. 
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Figure 1. RXTE- ASM and Swift-B AT light curves of 4U 
1626-67, binned with 30 d. The 15-50 keV Swift-BAT light curve 
is shown with a thick hne over the 1.5-12 kcV B.XTE- ASM Ught 
curve (thin line). The inset figure shows the expanded view near 
the onset of the flux enhancement which occured around MJD 
54500 (thick line is the Swift-BAT light curve in the same units 
as in the main panel). 



2 OBSERVATIONS AND ANALYSIS 

4U 1626—67 is a medium intensity persistent X-ray source 
among the few hundred bright X-ray sources regularly mon- 
itored by the AU Sky Monitor (ASM) on board Rossi X- 
ray Timing Explorer (RXTE). Figure 1 (thin line) shows 
the long term 1.5-12 keV ASM light curve of 4U 1626-67, 
binned with 30 d. The data used here covered the time range 
between MJD 50088 to MJD 54909. The X-ray intensity was 
gradually decreasing until MJD 54500, when an increase in 
the X-ray flux was noticed. The inset in Figure 1 shows 
the expanded view near the onset of the flux enhancement. 
There seems to be two episodes of X-ray flux enhancement, 
both on timescales of ~100 days. 

4U 1626—67 was also regularly monitored by the Burst 
Alert Telescope (BAT; Barthelmy et al. 2005) on board 
the Swift observatory (Gehrels et al. 2004). The long term 
15-50 keV Swift-BAT light curve, binned with 30 d is also 
shown in Figure 1 (thick line) alongwith the ASM light 
curve. The observations covered the time range from MJD 
53414 to MJD 54920. A sudden increase in the X-ray flux 
was first reported by Krimm et al. (2008). The inset in Fig- 
ure 1 (thick line) shows the expanded view (in the same 
units) near the onset of the flux enhancement. The flux 
change on timescales of ~100 days is also clear in the Swift- 
BAT light curve. 

We have also analysed data obtained from the Propor- 
tional Counter Array (PCA) and the High Energy X-ray 
Timing Explorer (HEXTE) on board the RXTE. The PCA 
consists of five xenon/methane proportional counter units 
(PCUs) and is sensitive in the energy range of 2—60 keV 
with an effective area of 1300—6500 cm^, depending on the 
number of PCUs ON (Jahoda et al. 1996). The HEXTE op- 
erates in the energy range 15—250 keV and consists of two 
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Table 1. Log of HXTE olwor vat ions ariah'zcd in tiio present work 



Yeaj: Observation ID Number of Total Exposure 
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1996 
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153 
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10 




P20146 
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1997 


P20146 


12 
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1998 


P30058 


2 


40 




P30060 


7 


53 


2008 


P93431 


2 


7 


2009 


P94423 


3 


9 



clusters of phoswich scintillation detectors, each having a 
collecting area of 800 cm'^ (Grubcr ct al. 1996). 

Table 1 gives a log of the RXTE observations used for 
the present analysis. We have studied the temporal and spec- 
tral variations in 4U 1626—67, using data collected in the 
Standard-1 and the Standard-2 mode of the RXTE-VCK, 
respectively. The archive mode data of the HEXTE was used 
for the spectral analysis. 



2.1 Timing Analysis 

We searched for pulsations using data obtained in the Stan- 
dard 1 mode of the RXTE-VCK. The background counts 
were simulated using the ftool - runpcabackest and sub- 
tracted from the source light curve. The photon arrival times 
wore converted to the solar system barycenter using the 
FTOOL - FAXBARY. We Searched for the spin period of the 
neutron star using the efsearch tool of the HEAsoft analy- 
sis package, ftools ver 6.5.1. This tool folds the light curve 
with a large number of trial periods around an approximate 
period and determines the best period by the maximiza- 
tion technique. Using this method, we determined a pulse 
period of 7.67941(1) s at MJD 54530.4, 7.67945(7) s at MJD 
54538.1 and 7.67848(2) at MJD 54984.5 (Jain & Paul 2009). 
This gives a spin-up rate of 4.02(5) x 10~^^ Hz s~^. Simi- 
lar behaviour based on the Formi/GBM data has also been 
reported by Camero-Arranz et al. (2009). 

Figure 2 shows the normalized pulse profile of 4U 
1626—67 binned into 64 phasebins. The top panel shows 
the pulse profile generated from data obtained from RXTE- 
PCA observations made in 1996 (ObsID 10101-01-01-00) 
and folded with a period of 7.66735 s (Jain et al. 2007). 
The bottom panel shows the normalized pulse profile of 2008 
data (ObsID 93431-01-01-00) folded with a period of 7.67941 
s. Although the two profiles resemble the characteristic bi- 
horned profile of the source (Lcvinc ct al. 1988), the pulse 
fraction is significantly different in the two cases. Moreover, 
the peaks in the 2008 pulse profile are significantly sharper 
with a higher amplitude variation than those in the 1996 
pulse profile. 

We searched for QPOs in the 2—60 keV energy band, 
using the data collected in the Standard 1 mode of PCA, 
which has a time resolution of 0.125 s. We created the power 
density spectrum (PDS) using the ftool - powspec for 
all the sufficiently long RXTE observations (See Table 1). 
The PDS were normalized such that their integral gives the 
squared rms fractional variability and the white noise level 
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Figure 2. The 2-60 keV normalized pulse profile of 4U 1626-67 
generated from observations made with RXTE-VCA during 1996 
(ObsID 10101-01-01-00) and 2008 (ObsID 93431-01-01-00). The 
background subtracted light curves were folded into 64 phasebins 
with a period of 7.66735 s and 7.67941 s, respectively. 

was subtracted. Figure 3 shows some representative power 
spectra generated from some of the observations mentioned 
in Table 1. The dates of observation are mentioned in each 
panel above their respective observation IDs. 

The narrow main peak at around 0.130 Hz corresponds 
to the spin frequency of the neutron star. Multiple harmon- 
ics are also seen in the PDS of the source. In addition to 
the main peak, a QPO feature is also seen at ^^48 niHz. 
However, the QPO feature is not present in the PDS gener- 
ated from observations made in 2008 and 2009 (Figure 3 (9) , 
(10), (11), (12)). In all the observations prior to 2008, the 
power spectrum consists of a broken power-law. As opposed 
to this, there is a remarkable difference in the power spectra 
generated from 2008 and 2009 observations. 

Figure 4 shows the power spectrum generated from the 
2008 observations (ObsID 93431-01-01-00; Figure 3-panel 9) 
overlaid with a typical power spectrum from a 1996 observa- 
tion (ObsID 10101-01-01-00; Figure 3-panel 1). In the 2008 
power spectrum, absence of a QPO at ~48 mHz is clearly 
seen. One more feature that stands out is the substantial 
broadening of the wings of the pulse peak in the power spec- 
trum during the 2008 observation. The r.m.s variability in 
the low frequency noise component of the power spectrum 
was 14.6%, whereas, the r.m.s variability in the broad wings 
of the fundamental pulse peak and the fundamental pulse 
peak itself was 6.9% and 5.2%. This broad feature at ^ 130 
mHz was never detected before and hence provides the first 
instance of a strong coupling between the aperiodic and pe- 
riodic X-ray variability in the source. 
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Figure 3. Some representative X-ray power spectra of 4U 1626—67 generated from tlie RXTE-PCA observations. The date of observation 
is mentioned in each panel above the respective observation ID. 



2.2 Spectral Analysis 

The pulse-phase averaged spectra of 4U 1626—67 were gen- 
erated from the data collected in the Standard 2 mode of the 
RXTE-PCA and the archive mode of the RXTE-UEXTE. 
The background counts for the PCA data were simulated us- 
ing the FTOOL - RUNPCABACKEST and subtracted from the 
source spectrum. In case of the HEXTE data, the source and 
the background spectra were generated using the ftool - 
HXTBACK. These were then corrected for the detector dead 
time using the ftool - hxtdead. In order to bring out the 
significant evolution in the X-ray spectrum, wc have gener- 
ated the spectrum from one of the 1996 observation (Ob- 
sID 10101-01-01-00) and the 2008 observation (ObsID 
93431-01-01-00). 



Figure 5 shows the 3—20 kcV and 15—50 kcV energy 
spectrum of the source from the 1996 RXTE-PCA and 
RXTE-BEXTE observations. A simultaneous fit of the two 
spectrum consists of an intrinsic absorption, a powerlaw, a 
gaussian, a high energy cutoff and cyclotron resonance fea- 
ture (CRF) . The bottom panel of Figure 5 shows the residual 
of the fit. The best fit spectral parameters (with 1% system- 
atic errors) are given in Table 2. The spectrum is similar 
to the broadband continuum spectrum obtained from the 
Beppo-SAX. observations made after the torque reversal in 
1990 (Owens et al. 1997; Orlandini ot al. 1998). The Beppo- 
SAX energy spectrum was well fit by a model consisting of a 
low energy absorption, a blackbody, a powerlaw and a high 
energy cutoff. However, the best fit blackbody temperature 
was ^0.3 keV which is much lower than the sensitivity of 
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Figure 4. X-ray power spectra of 4U 1626-67 from the RXTE 
observations of 1996 and 2008. The upper, dotted curve is the X- 
ray spectrum of 1996 observation (ObsID 10101-01-01-00; Figure 
3-panel 1) and the lower, solid curve is from 2008 observation 
(ObsID 93431-01-01-00; Figure 3-panel 9). 




Energy (keV) 

Figure 5. The X-ray spectrum of the 1996 RXTE-^CK and 
RXTE-HEXTE observations of 4U 1626-67. The model compo- 
nents are intrinsic absorption, powerlaw, a gaussian and a high 
energy cutoff. The bottom panel shows the residual of the fit. 

the RXTE-VCA detector. We have obtained a high energy 
cutoff of 6.45 keV which is also quite small as compared to 
~20 keV obtained from the Seppo-SAX observations. 

The 3-20 keV RXTE-VGA and the 15-50 keV RXTE- 
HEXTE energy spectrum obtained from the 2008 observa- 
tions are shown in Figure 6. The spectra were simultaneously 
fitted with a model consisting of a powerlaw, a blackbody, a 
gaussian, high energy cutoff and cyclotron resonance feature 
(CRF). In Table 2, we have given the best fit spectral param- 
eters (with 1 % systematic errors). The observations made 
with HEAO-1 (Pravdo ct al. 1979) and X-ray astronomy 
satellite Tenma (Kii et al. 1986) during the earlier spin-up 
phase of 4U 1626—67 also showed similar energy spectrum. 

Figure 7 shows the evolution of the 3—50 keV X-ray 




5 10 20 

Energy (keV) 



Figure 6. Energy spectrum of 4U 1626—67 obtained from the 

2008 RXTE-PCA observations. The model consists of a black- 
body, a powerlaw and a gaussian component. 



Table 2. Best fit spectral parameters of 4U 1626—67, obtained 

from the RXTE-PCA and BXTE-llEXTE observations. 



Parainclci' \"ahi(* 





1996 


2008 


Nh (10^2 cm-2) 


0.65 


0.82 


Photon index (F) 


0.59(4) 


0.74(7) 


norm (photons keV"'^ cm"^ ^~^) 


0.0079(6) 


0.012(2) 


Gaussian line (kcV) 


6.56(40) 


6.33(39) 


Gaussian width (keV) 


1.15(21) 


1.43(38) 


Gaussian norm (photons cm^'^ 


0.00054(32) 


0.0019(18) 


Blackbody, kT (keV) 




0.66(12) 


Blackbody, norm (km-^ kpc"'^) 




44(7) 


cutoff E (keV) 


6.45775(95) 


17.7(8) 


foldE (keV) 


32(5) 


26(12) 


Energy(cflF) (keV) 


39.5(1.8) 


40.3(1.2) 


Width(CHF) (keV) 


13.8(3.0) 


4.8(4.5) 


Depth(CHF) 
Reduced (d.o.f.) 


2.02(28) 


2.9(1.7) 


0.59 (57) 


0.61 (43) 



spectrum from 1996 observation to the 2008 observation. 
In the figure, the best fit spectrum of the 1996 observation 
(dotted curve) is shown with the spectrum generated from 
the 2008 observations. The bottom panel in the figure gives 
the ratio between the two spectra. It is clear that there is an 
excess soft-photon emission in the spectrum generated from 
the 2008 observations. The powerlaw photon index has also 
increased during the recent observation. 



3 DISCUSSION 

4U 1626—67 is a unique system in which the X-ray flux has 
been gradually decreasing over the past ~30 years. Since 
the discovery of X-ray pulsations in 1977, the source has 
shown a smooth spin evolution history. Initially, it was found 
to be spinning up with a characteristic timescale of 5000 
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Figure 7. The best fit spectrum of 4U 1626—67, generated from 
the 1996 (dotted line) and 2008 (solid line) RXTE observations. 
The bottom panel shows the ratio of the two spectra. 



yr, but after an abrupt torque reversal, the neutron star 
started spinning down at about the same rate (Chakrabarty 
et al. 1997; Krauss et al. 2007). The recent flux enhancement 
has opened up many questions regarding the stability of the 
accretion disk. We have found a sudden increase in the X-ray 
flux of 4U 1626—67, which could have also caused significant 
changes in the geometry of the accretion column. 

The power spectra of X-ray pulsars consist of narrow 
peaks originating from the periodic signal and its harmonic 
lines which are often accompanied with aperiodic features 
like broad bumps, wiggles and steep low frequency excess 
(Wijnands & van der Klis 1999). The periodic variabilites 
arise due to the rotation induced motion of the accretion 
column, through which the infalling matter is funneled to- 
wards the magnetic poles of the neutron star. Whereas, any 
instability in the emissivity of the accretion column, can give 
rise to an aperiodic variability (Burderi et al. 1997). If the 
inhomogeneities occur far away from the surface of the neu- 
tron star, then the periodic and aperiodic variations occur 
independent of each other (Angelini et al. 1989). But if the 
aperiodic variability occurs near the accretion column, then 
they are modulated at the pulsar spin frequency (Makishima 
et al. 1988). Lazzati & Stella (1997) found a significant cou- 
pling between the aperiodic and periodic variabilities in Vela 
X-1 and 4U 1145—62. A significant coupling was also re- 
ported in the accretion powered millisecond binary pulsar 
SAX J1808.4-3658 (Menna et al. 2003), Her X-1 (Moon 
& Eikenberry 2001) and Cen X-3 (Raichur & Paul 2008). 
The coupling between the aperiodic and periodic variations 
is expected to alter the shape of the power spectrum be- 
cause of the convolution of the two signals (Burderi, Robba 
& Cusumano 1993). The coupling results in a broadening of 
the base of the power spectrum peaks, as if, the peats are 
riding on a QPO. 

X-ray spectrum and rapid X-ray variability of accreting 
compact objects have a common origin and are therefore ex- 
pected to be correlated. In the case of 4U 1626—67, Yi & 
Vishniac (1999) had found that the torque reversal in 1990 



was accompanied by spectral transition and a change in the 
luminosity. During the spin-up phase of 4U 1626—67, the 
phase averaged spectrum was well fit by a blackbody tem- 
perature kT of ~0.6 and a powerlaw photon index of ~1 
(Pravdo et al. 1979; Kii et al. 1986). Vaughan & Kitamoto 
(1997) found that the spin transition had a significant effect 
on the 2—10 keV energy spectrum. The powerlaw photon in- 
dex changed from 1.6 (Pravdo et al. 1979) to ~0.6 (Owens et 
al. 1997). The blackbody temperature also decreased from 
~0.6 keV to ^0.3 keV. We have found similar spectral tran- 
sition in the present high flux state of the source. The spec- 
tral parameters derived during the present high flux state 
arc similar to those determined during the earlier spin-up 
phase. 
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